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used  for  any  purpose  other  than  in  connection  with  a  definitely  related 
Government  procurement  operation,  the  United  States  Government  there¬ 
by  incurs  no  responsibility  nor  any  obligation  whatsoever;  and  the  fact 
that  the  Government  may  have  formulated,  furnished,  or  in  any  way  sup¬ 
plied  the  said  drawings,  specifications,  or  other  data,  is  not  to  be  re¬ 
garded  by  implication  or  otherwise  as  in  any  manner  licensing  the  hold¬ 
er  or  any  other  person  or  corporation,  or  conveying  any  rights  or  per¬ 
mission  to  manufacture,  use,  or  sell  any  patented  invention  that  may  in 
any  way  be  related  thereto. 

The  informationfurnishedherewith  is  made  available  for  study 
upon  the  understanding  that  the  Government’s  proprietary  interests  in 
and  relating  thereto  shall  not  be  impaired.  It  is  desired  that  the  Jidge 
Advocate  (WCJ),  Wright  Air  Development  Center,  Wright -Patterson  Air 
Force  Base,  Ohio,  be  promptly  notified  of  any  apparent  conflict  between 
the  Government’s  proprietary  interests  and  those  of  others. 
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FOREWORD 

This  report  was  prepared  "by  the  Mapping  and  Charting  Research 
Laboratory  of  the  Ohio  State  University  Research  Foundation,  under 
USAF  Contract  Ho,  AF  18(6oo)90,  The  contract  was  administered  under 
the  direction  of  the  Mapping  and  Charting  Branch,  Photographic 
Reconnaissance  Laboratory,  Air  Research  and  Development  Command, 
Wright-Patterson  Air  Force  Base,  Dayton,  Ohio  with  Mr,  D.  L,  Radcliffe, 
Chief  of  the  Mapping  and  Charting  Branch,  as  Project  Engineer, 

Research  and  Development  Order  Ho,  R683-44,  "Charting,  Aeron¬ 
autical,  Photogrammetry  and  Geodesy,"  and  R683-58,  "Aeronautical 
Charting  Systems,"  are  applicable  to  this  report. 

This  report  is  identified  by  the  Ohio  State  Research  Foundation 
as  03URF  Technical  Paper  156,  project  485. 
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ABSTRACT 


An  analytical  method  of  determining  the  exterior  orientation 
of  a  near-vertical  photograph  is  described.  The  theory,  computing 
procedure,  test  results,  and  sample  computations  using  fictitious 
photographs  are  presented.  It  is  concluded  that  the  method  is  rapid 
and  meets  practical  accuracy  needs  for  photographs  containing  tilts 
up  to  three  degrees. 

PUBLICATION  REVIEW 

This  report  has  been  reviewed  and  is  approved. 

FOR  THE  COMMANDING  GENERAL: 

l  GORDON  A.  BLAKE 

Brigadier  General,  USAF 

Chief,  Weapons  Components  Division 
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EXTERIOR  ORIENTATION  OF  NEAR-VERTICAL  AERIAL  PHOTOGRAPHS  - 
COMPUTATIONAL  PROCEDURE  AND  ILLUSTRATED  EXAMPLES 

INTRODUCTION 

The  analytical  method  of  determining  the  elements  of  exterior 
orientation  of  an  aerial  photograph  presented  here  is  unique  in  that  a 
complete  determination  of  tilt  is  made  without  computing  the  space 
position  of  the  photographic  exposure  station.  Eliminating  the  necessi¬ 
ty  for  the  space  position  determination  greatly  simplifies  the  computing 
procedure. 

The  theory  for  this  method  of  determining  tilt  was  developed  by 
Mr.  W.  0.  Byrd  and  presented  as  part  of  Technical  Paper  No.  li*2, 

"Some  Elementary  Aspects  of  Computational  Problems  of  Photograametry, " 
prepared  at  this  Laboratory.  In  that  paper  he  developed  several  methods 
of  determining  the  elements  of  exterior  orientation  and  exposure  station 
position,  or  combinations  of  both.  The  test  results  and  computational 
procedures  of  these  methods  will  be  subjects  for  subsequent  papers.  In 
this  paper  the  theory  presented  by  Mr.  Byrd  is  restated. 

The  computational  procedure  for  the  method  presented  here  is 
shown  as  a  computing  form.  Sample  computations,  fictitious  photographs 
for  test  purposes,  and  a  guide  form  indicating  all  computing  operations 
are  included. 
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SECTION  I 


THEORY 

1.1  A  restatement  of  the  tneory  pertaining  to  this  method  of  tilt 
analysis  is  required  for  a  thorough  study  of  the  computational  pro¬ 
cedure  and  final  results.  The  evolution  of  the  theory  and  the  equations, 
as  presented  here,  have  been  taken  directly  from  Technical  Paper  No.  11*2 
by  Mr.  W.  0.  Etyrd  of  this  Laboratory. 

1.2  The  direction  cosines  relating  the  photographic  coordinate  system 
with  the  ground  coordinate  system  are  tabulated  as  follows: 

cos  x  cos  y  cos  z 

X  kx 


These  nine  elements  comprise  the  basic  concept  governing  the  orientation 
between  the  two  coordinate  systems. 

1.3  Consider  now  a  directed  line,  PL,  that  intersects  the  ground  plane 

at  point  P(Xp,  Yp,  Zp),  the  photo  plane  at  p(x^,  y^),  and  the  exposure 

station  at  L(X  ,  Y  ,  Z  ).  Then  the  direction  cosines  of  this  line  with 
L  L  L 

respect  to  each  coordinate  system  are  expressed  as  follows: 

with  respect  to  the  Ground  System 
cos  X  cos  Y  cos  Z 

PL  Px  P2  P3 

and 

with  respect  to  the  Photographic  System 
cos  x  cos  y  cos  z 

PL  Px  Pa  P3 
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where 


h’h> 

pi - =•“  “  “iPi  +  %Pa  +  *iP3 

Y.  -Ip 

p  -  .!*_  f  -  nijP^  +  n2p2  +  k2p3  (1-1) 

PL 

z  —  z 

P3  -  J^£  -  m3pi  *  “A  ♦  k3p3 

* 

and 

p  -  -  -  n^Pj.  +  m^  +  m3P3 

pL 

Pa  -  -  »ipi  *  napa  +  n3p3  (I“2^ 

Pa  -  *jz’  ^  *  k*P»  *  kaPa  • 

1.1*  From  the  above  relationships,  equations  may  be  written  for  the 
determination  of  either  the  photographic  coordinates  or  the  ground 
coordinates  in  terms  of  the  other  coordinate  system  and  the  nine  direction 
cosine  elements.  By  ratioing,  substituting,  and  reducing  the  above 
relationships,  the  photographic  coordinates  are  expressed  thus, 


Xg  _  *1%  *  *a*P  *  a3zp  ♦  "4 

f  kjXp  ♦  kaYp  ♦  k3Zp  ♦  k4 

Yg  _  r^Xp  +  HgYp  +  n3Zp  ♦  n4 

f  kjXp  ♦  kaYp  +  k3Zp  +  k4 


(1-3) 
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where 


®4  *  -  +  m*YL  +  ®3Zl) 

»4  "  ~  <Vl  *  n2YL  +  n3ZL) 

*4  -  -  <kiXL  +  k2YL  +  k3ZL)  * 

Rewriting  the  above  equations,  we  have 

k, 

10j.Xp  ♦  m2^p  +  ®3Zp  ♦  BI4  +  (Xp  Xp) 

+  T  (xP  V  +  T  (xpzP>  +  T  xp  "  0  (**> 

and 

k. 

niXP  +  n2YP  +  n3ZP  +  n4  *  T"  (ypXp) 

♦  t  (ypYp) +  t  (ypV  *TxPm  0  •  (I_5) 

1.5  In  order  to  eliminate  m4  and  n4  ,  substitute  into  I-U  and  1-5 
the  coordinates  of  two  points.  Let  the  ground  coordinates  of  the 
points  be  A(X^,  ^A»  ZA)  and  B(Xg,  Yg,  Zg),  and  the  photographic 
coordinates  are  then  aCxj,  ya)  and  b(x^,  y^) .  Subtracting  the  two 
equations  will  form  the  following: 

(Xp  -  XA)  +  m2(YB”YA^  +  %(ZB  ”  ZA> 

+  T  (xbXB“xaXA)  +  T  (Vb-Va) 

K 

+  f  (xbZB“  xaZA^ 

k4 

+  T(Vxa)*°  (1-6) 
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and 


°i(*B  ”  YA^  +  +  n3^ZB-^A^ 

♦  £  1^%-^*)  *  T  I'b1!-*) 

k  k 

♦  f  (?bV  *aV  *  -f  (/b-i’a)  *  0  •  U-7) 

j£ 

Eliminating  by  multiplying  (1-6)  by  (y^-y^  and  (1-7) 
by  (x^-Xj^,  then  subtracting  (1-7)  from  (1-6),  we  have 

%(yb-ya)(xB-xA> +  m2^b-yaHYB^  V 

♦  “»3(yb  -yaH^B"  ZA)  -  nj,  (xb -  xa)  ( X0  «-  XA) 

“na(xbvXa)<YB-IA)  -  “a^^W 

♦  yC*eyb  "Vj)  (XB"Xa) 

+  /  <Vb-V«>  <yb~ya> 

k 

♦"T  (xaXb“ybXa)  (zb~za)  -  0  *  (1-8) 


1.6  For  near- vertical  photography,  the  elements  of  exterior  orienta¬ 
tion  are  defined  as 


where 


(a-s)  ■  180  ♦  P 

V 


4 
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then,  the  nine  direction  cosines  relating  the  coordinate  systems  in 
terms  of  the  small  angles  p ,  t^,  And  ty  are  determined  and  tabulated 
thus. 

Near  Vertical 
(a  -  s)  -  180  ♦  p 


COS  X 

cos  y 

cos  z 

H 

X 

2  2 

+  ty  ♦  P  tat 

ft 

Y 

2 

♦1-4-4- 

2  2 

♦  tx  -  P  ty 

Z 

"  V 

“  tx 

♦i-iL-iL 
2  2 

The  complete  derivation  for  each  of  these  elements  to  terms  of  the 
second  order  is  shown  in  Technical  Paper  11*2,  pages  111-112. 
Substituting  the  near-vertical  orientation  elements  in  equation  1-8, 
grouping  the  first  and  second  order  terms,  and  letting 

XAB  "  XB  “  XA 


ZAB  *  h  “  ZA 
xab  ^  “  xa 

yat "  -  ya 


wab  "  xa  *b  -  *b  y* 
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gives  the  final  equation  as 


v  wab  YAB 

“  (xab  XAB  +  yab  YAB^  +  ^xab  ZAB  +  *  )  Sc 


-ab  XAB 


♦  (JiH-JS.  -  yab  Zab)  ty 

+  (yab  XAB  "  xab  YAB  +  f 


)♦  G 


AB 


(1-9) 


where 


*ab  YAB 


,%b  XAB  _  . 

°AB  "  +  (  f  )  t  Sc  "  (  f 


)  fi  ty 


*  <■**,  r*B  -  W  £ 


+  <*ab  r*B  * 


*ab  XAB.  ^x 
f  '  2 


^ab  XAB 

-  (yab  XAB  ♦  - T^)  2 


Sc  v 

+  (xab  XAB  _  yab  YAB'  2  • 

The  minimum  data  required  to  solve  for  p,  t^,  and  ty  is  three  control 
points.  Given  the  coordinates  of  these  points,  three  linear  non- 
homogeneous  equations  can  be  formed  from  1-9  by  letting  Gab  equal  zero 
in  the  first  solution.  If  a  second  solution  is  required,  the  constant 
terms  in  the  equations  can  be  corrected  by  using  the  first  solution 
values  for  p,  t^,  and  ty  and  computing  and  applying  Gag  as  illustrated 
in  the  appendix  (page  17). 

1.7  If  desired,  the  space  position  of  the  exposure  station  can  easily 
be  determined  by  the  following  sequence  of  equations: 
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* 


or 


,  .2  r  m3<-  V  *  n»(~  *  *>f 

L  P  r  m*  (  -  xp)  +  n*  ( -  yp)  +  f 

a3(-  ^p)  *  n3(-  yp>  *  k3f 

ZL  "  P  P  ma(-  Xp)  ♦  n2(-  yp)  ♦  ka  f  . 


(1-10) 


Using  the  computed  2^  the  other  coordinates  of  the  exposure  station 
are  then  determined  by 


%(-  xp)  ♦  yp) *  kx  f 

^  -  Xp  ♦  (Zl  -  Zp)  m3(-  Xp)  +  n3{-  yp) ?  k3  f 

Y  .  Y  WL-Z)  *»("  *»>>  *  V  1  f 

YL  P  ^  P  m3(-  Xp)  +  n3(-  yp)  ♦  k3  f 


(1-11) 


For  any  further  investigation  requiring  the  use  of  the  photographs,  it 
will  be  necessary  to  compute  the  nine  direction  cosine  elements  from 
the  determined  values  of  exterior  orientation.  Knowing  these  values, 
the  coordinates  of  the  space  position  can  readily  be  determined  from 
1-10  and  I- 11.  The  equations  for  evaluating  the  nine  cosines  are 
given  on  page  25  of  Technical  Paper  No.  11*2. 
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SECTION  II 


FICTITIOUS  PHOTOGRAPHS 

2.1  In  order  to  test  this  method  for  determining  tilt*  a  set  of  three 
overlapping  fictitious  photograpns  were  prepared.  The  overlapping 
condition  is  not  needed  in  this  case*  but  was  incorporated  for  use  in 
future  tests.  In  Figure  1  the  exterior  orientation  elements  and  space 
position  of  each  photograph  are  listed  with  a  diagram  showing  the 
approximate  position  of  the  computed  photo  points.  In  Table  I  the  ground 
coordinates  and  corresponding  computed  photographic  coordinates  of  each 
point  are  listed. 

Photo  1 

t  -  1*  00'  ^  «  15,000 

s  -  60*  00*  yl  -  15,000 

a  -  200*  00'  ZL  -  20,000 

Photo  2 

t  -  3*  00'  xL  -  25,000 

s  -  120*  00*  Yl  -  lli,  900 

a  -  260*  00'  ZL  -  20,100 

Photo  3 

t  -  5*  00*  -  35,000 

s  -  120*  00'  IL  -  15,200 

a  -  260*  00'  ZL  -  20,200 

Figure  1  -  Orientation  Elements  and  Exposure 
Station  Coordinates  in  Feet 
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Table  1  -  Photo  and  Ground  Coordinates 


2.2  A  focal  length  of  150  mm  was  used  for  all  fictitious  photo 
computations,  and  the  photo  coordinates  were  computed  by  the  method 
presented  in  Technical  Paper  No.  34,  "Fictitious  Photographs  for  aesearch," 
by  W.  0.  Byrd  of  this  Laboratory. 
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SECTION  III 


COMPUTING  FORM 

3.1  The  final  equation  (1-9)  on  which  this  method  is  based  can  be 
conveniently  adapted  to  a  computing  form  which  clearly  illustrates  the 
computing  procedures.  A  form  for  the  complete  tilt  determination  is 
shown  in  the  appendix  (page  17).  As  noted  in  the  form,  an  attempt  has 
been  made  to  place  and  group  all  values*  terms,  and  operations  in  the 
most  advantageous  position  for  the  computer's  use. 

3.2  All  but  one  of  the  computing  operations  are  shown  on  the  form. 

This  one  exception,  the  solution  of  the  three  non-homogeneous  linear 
equations,  can  be  accomplished  hy  several  different  methods.  In  this 
particular  case,  the  solution  was  performed  by  determinants  according 
to  Cramer's  rule.  By  this  method  the  second  solution  values,  as 
indicated  on  the  form,  can  be  conveniently  obtained  and  with  a  minimum 
number  of  computing  operations.  In  order  to  illustrate  the  method  used 
in  solving  the  three  operations,  a  guide  form  (on  page  18)  was  prepared. 
Since  the  computing  form  indicates  all  other  operations,  the  guide  form, 
in  most  cases,  will  be  used  by  the  computer  only  to  become  familiar  with 
the  equation  solution.  However,  all  operations  are  shown  on  the  guide 
form. 

3.3  According  to  the  assumption  (u  -  s)  »  180  ♦  £,  angle  {I  must  be 
smallj  therefore,  it  may  be  necessary  to  rotate  the  photo  coordinates 
through  an  aigle  3.  The  angle  3  can  be  determined  computationally  or 
graphically,  whichever  is  preferred,  and  the  rotation  of  coordinates  is 
then  performed  as  indicated  on  the  form.  If  3  is  determined  graphically, 
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the  plot  can  be  used  to  check  approximately  the  magnitude  and  sign  of 
the  rotated  photo  coordinates.  To  determine  0  graphically,  a  plot  of 
the  ground  points,  at  some  convenient  scale,  is  prepared  on  vallum  or 
tracing  paper.  On  the  photo,  radial  lines  from  the  principal  point  are 
drawn  through  the  image  points.  The  vellim  plot  of  the  ground  points 
is  superimposed  on  the  radial  lines  of  the  photo  images*  0  is  the  angle 
between  the  photo  and  ground  x-axes. 

3.k  For  convenience  of  computation,  the  factor  of  10  3  was  used. 

Since  each  term  of  the  basic  equation  (1-9)  contains  a  ground  coordinate 
value,  the  introduction  of  this  factor  has  no  effect.  Also,  in  the 
formation  of  the  three  equations,  the  factor  merely  displaces  the 
decimal  point  in  all  terms. 
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SECTION  IV 


TEST  RESULTS 


it.  1  Computations  were  made  for  each  of  the  three  fictitious  photo¬ 

graphs,  with  angle  p  the  same  magnitude  as  the  tilt*  also,  computations 
were  made  for  Photo  3a,  with  angle  p changed  to  one  degree.  A  complete 
set  of  computations  are  included  in  the  appendix.  The  results  of  the 
computations  were  tabulated  as  follows: 


Photo  1 


Actual  Values 

Computed  Values 

First 

Solution 

Error 

Second 

Solution* 

Error 

ym 

l 

8 

0*  59.9' 

-  o.l* 

t  -  1*  00* 

0*  59.5' 

-  0.5' 

s  -  60*  00* 

59*  39.1' 

-  20.9' 

a  •  200*  00' 

199*  39.2' 

-  20.3' 

0  -  39*  00* 

Photo  2 

Actual  Values 

Computed  Values 

First 

Solution 

Error 

Second 

Solution* 

Error 

P  -  y  00* 

3*  00.2' 

♦  0.2* 

2*  59.8' 

-  0.2' 

t  -  3*  00' 

2fc  57.6* 

• 

CM 

1 

2*  59.U' 

-  0.6' 

s  -  120*  00' 

119*  1*3.8' 

-  16. 2» 

119*  59.2* 

-  0.8' 

u  -  260*  00* 

259*  1*3.6* 

-  16.1*' 

259*  59.1*' 

-  0.6' 

e.«  37*  oo' 

♦Second  order  term  solution 
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Photo  3 


Actual  Values  Computed  Values 


First 

Second 

Solution 

Error 

Solution* 

Error 

fi  -  5* 

00* 

1*#  58.1* 

1.9* 

1**  59.U* 

0.6* 

t  -  5* 

00* 

5*  29.1* 

+  29.1' 

k*  53.1*' 

6.6* 

a  -  120* 

00* 

113*  35.1* 

♦  6*  21*. 9* 

117*  1*0.0* 

-  2*  20.0* 

a  -  260* 

00* 

253*  37.0* 

-  6*  23.0* 

257*  l*o.6» 

-  2*  19.1*' 

0  -  37- 

00* 

Photo  3a 

Actual  Values 

Computed  Values 

Firet 

Solution 

Error 

Second 

Solution* 

Error 

P  * 

1* 

00* 

0*  57.0' 

3.0* 

1*  00.2* 

♦ 

0.2* 

t  » 

5* 

00' 

5*  29.1* 

♦  29.1' 

1*#  53.1*' 

- 

6.6* 

s  - 

120* 

00* 

117*  1*1.6* 

-  2*  18.2* 

120*  27.1*' 

♦ 

27.1*' 

a  ■ 

260* 

00* 

257*  1*1*.  3' 

-  2*  15.2* 

260*  27.2* 

♦ 

27.2* 

9  - 

39. 

00* 

■••Second  order  tern  solution 
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In  Table  II  only  the  errors  of  the  final  results  are  tabulated 
for  comparison. 


Photo 

Actual  Values 

Errors 

-  Computed  Values 

t 

s 

P 

t 

s 

1 

1*  00* 

1*  00' 

60*  00' 

-  0.1* 

-  o.5» 

20.9' 

2 

3*  00* 

3*  00' 

120*  00' 

-  0.2* 

-  0.6' 

0.8* 

3 

5*  oo* 

5*  00' 

120*  00' 

-  0.6' 

-  6.6« 

-  2*  20.0* 

3a 

1*  oo* 

5*  oo* 

120*  00' 

♦  0,2* 

-  6.6* 

♦  27.ii* 

Table  II 

-  Errors 

in  Final  Results 

Since  a  (a  -  180  ♦  0  ♦  s)  is  defined  in  terms  of  p  and  s, 
the  magnitude  of  .error  in  a  will  be  the  sum  of  errors  in  p  and  s. 
Hence,  in  Table  II  a  was  emitted. 


VADC  TR  52-95 


15 


SECTION  V 


CONCLUSIONS 

5.1  The  development  of  the  general  basic  theory,  which  is  presented 
in  Technical  Paper  No.  li±2,  should  be  emphasized  because  of  its 
significant  importance  to  research  problems  in  which  photographs  are 
used.  In  this  particular  case,  the  theory  has  given  rise  to  a  practical 
and  rapid  analytical  method  of  determining  the  photographic  elements  of 
exterior  orientation. 

5.2'  For  photographs  with  tilts  up  to  three  degrees,  as  indicated  by 
the  results  snown  in  Table  II,  the  method  meets  practical  accuracy 
needs  for  a  determination  of  exterior  orientation  elements.  Since  the 
majority  of  aerial  photography  flown  todsy  contains  less  than  3*  00' 
and  usually  closer  to  1*  00'  tilts,  this  method  is  applicable  and  will 
prove  to  be  a  useful  tool  in  pho tog ranme trie  research  and  mapping 
problems. 

5.3  The  time  required  for  performing  a  complete  solution,  including 
second  order  terms,  is  about  two  hours.  However,  due  to  the  small 
tilts  in  a  majority  of  the  photographs,  the  second  solution  will  be 
unwarranted  and  the  computing  time  will  be  reduced  considerably. 

5.U  The  results  of  Photo  3a  show  that  an  improvement  of  U*  00'  in 
{3  angle  has  little  effect  on  the  resulting  tilt  angle,  but  does  increase 
the  accuracy  of  the  swing  angle. 

5.5  The  method  as  presented  does  not  provide  a  check  on  the  cotnpur- 
tations.  However,  the  completed  computations  can  be  checked  by  sub¬ 
stituting  p,  tj,  and  ty  in  the  original  equations. 
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